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A Fermi Liquid theory is developed for the persistent current past a side coupled quantum dot 
yielding analytical predictions for the behavior of the first two harmonics of the persistent current 
as a function of applied magnetic flux. The quantum dot is assumed weakly coupled to a ring 
of non-interacting electrons and thus appropriately described as a Kondo impurity. The theory is 
valid at weak Kondo couplings in the regime where the system size, L, is much larger than the 
size of the Kondo screening cloud, £k. The predictions of the Fermi Liquid theory are compared to 
exact diagonalization results for the persistent current that lend support to the existence of a regime 
correctly described by this theory. The finite temperature conductance, at T <C Tk is also calculated 
using Fermi liquid theory allowing the definition of a "Wilson ratio" relating the conductance and 
the persistent current. 

PACS numbers: 72.10.Fk, 73.23.Ra, 72.15.Qm, 73.23.-b 



I. INTRODUCTION 

The recent experimental observation of the Kondo 
effect and related physical phenomena in quantum 
dotSfi^i^ quantum corrals 5 molecular electronics de- 
vices^ and carbon nanotubes^ as well as related nano- 
structures has led to significant renewed interest in this 
fundamental effect. For a review of recent progress see 
Ref. 1. The Kondo effect in ordinary metals is usu- 
ally associated with an increase in the resistance at low 
temperatures but experiments performed on these nano- 
structuresi^££ have shown that the Kondo effect can 
also lead to an increase in the conductance that in some 
settings can reach the unitary limit (G = 2e 2 /h), com- 
pletely overcoming the Coulomb blockade. The funda- 
mental and technological interest in the complete un- 
derstanding of this phenomenon is therefore consider- 
able. The experiments are typically performed on semi- 
conductor quantum dots where the electron occupation 
number on the dot is controllable by a gate voltage. 
When the occupation number is odd the quantum dot 
can act as a spin s = 1/2 impurity screened by the elec- 
trons in the leads. The screening of the impurity spin by 
the conduction electrons is, in some circumstances, asso- 
ciated with the formation of a "screening cloud" of size 
S.K = vf/Tk surrounding the impurity.— dU&liiJdilSJIi 
Here vf is the Fermi velocity and Tk the Kondo temper- 
ature. 

We note that the picture of a screening cloud of size 
vf/£,k is valid in a one dimensional system and also in 
an infinite two or three dimensional system where only 
one channel (such as the s-wave) couples most strongly 
to the Kondo impurity. However the behavior may be 
rather different in a disordered two or three dimensional 
box in which all states couple strongly to the Kondo im- 
purity. For the 3-dimensional model considered in [17j . 
for instance, the effective screening cloud is much smaller, 



oc T K 1/3 . See also [Hj]. 

The experimental observation of the screening cloud 
has proven elusive but recent experiments^ have sug- 
gested that it might be observable in finite-size prop- 
erties of the conductance when the dot is attached to a 
single quantum wire that forms an Aharonov-Bohm ring 
of length L. A natural extension of these experiments^ is 
to consider not only the conductance but also finite-size 
effects observable in the persistent current generated in 
the presence of a magnetic flux, $. 

The persistent current in quantum dot systems oper- 
ated in the Kondo limit has been the topic of a series 
of recent theoretical stn & es iML2L&&&&&£L22&aB. 
Usually, two fundamentally different geometries are 
considered'^ the embedded quantum dot (EQD) and 
the side-coupled quantum dot (SCQD) and it has been 
clearly established that the persistent current has a 
strong dependence on the parity of the total number of 
electrons, N, in the system. Here we shall exclusively 
be concerned with the SCQD. Suppressing electron spin 
indices, the Hamiltonian is in this case: 
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(c L =Cq) (1.1) 



where Hk = JkS ■ c' %Co, and Cj a is the electron anni- 
hilation operator at site j for spin a and S a are S=l/2 
spin operators. In the following we set t = 1. We always 
include the impurity spin as an electron in our definition 
of N but do not count the impurity site in our definition 
of L. Hence, at half-filling we have N = L + 1. The 
dimensionless magnetic flux, a = e<fr/c is introduced by 
appropriately changing the phases of the hopping terms 
at the corresponding sites. (We set h = 1.) The persis- 
tent current is given by j — —edEo/da where Eq is the 
ground-state energy. Despite the simplicity of this ge- 
ometry, a complete understanding of the behavior of the 
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persistent current both at half-filling and off half-filling 
has proved exceedingly difficult to achieve. In previous 
by one of us and P. Simon it was predicted 



25.2(i 



that, jL/(evp) — > when £,k/L — > for both parities of 
N in contrast to the results of Refis 123,24. However, sub- 
sequent work^ - confirmed this prediction at half-filling 
but showed both analytical and numerical evidence for 
additional terms off half-filling for N even. 

In the present paper we develop a Fermi liquid the- 
ory (FLT) for the persistent current past a SCQD, valid 
for I <C £,k <C L. The theory provides independent 
confirmation of the above mentioned prediction at half- 
filling 25,26 as well as additional analytical evidence for 
the additional terms off half-filling for N even. The the- 
ory yields detailed predictions about the finite-size in- 
dependence of the first two harmonics of the persistent 
current, using a single parameter, 1/Tr-. Section |TT] out- 
lines the FLT in considerable detail while section IIII Al 
presents our numerical results at half-filling and sec- 
tion IIII 51 our numerical results off half-filling. In Sec. IIVI 
we calculate the finite temperature conductance past the 
quantum dot, at T <C Tk ■ This enables us to calculate a 
generalized "Wilson ratio" relating the persistent current 
for a finite ring at T — with the conductance at finite T 
for the same quantum dot side-coupled to infinite leads. 
This could be useful in an experimental context since a 
prior infinite leads measurement of the low T conduc- 
tance would fix the precise numerical coefficient in the 
finite size dependence of the persistent current at T = 0. 



II. FERMI LIQUID THEORY FOR THE 
PERSISTENT CURRENT 



Kondo physics at low energy scales and long length 
scales, T <C Tk, L 3> £k, can be described by an ef- 
fective Hamiltonian not containing the impurity spin, 
since it is screened by the conduction electrons and the 
energy scale associated with breaking this spin singlet 
is of 0{Tk)- The low energy quasi-particles propogate 
freely except that they obey a modified boundary con- 
dition which can be thought of as arising because they 
must stay in wave-functions orthogonal to that of the 
screening electron. This modified boundary condition 
is equivalent to a 7r/2 phase shift for even-channel elec- 
trons (or for s-wave electrons in the three dimensional, 
spherically symmetric version of the model). In addition 
to the modified boundary condition an interaction term 
is generated in this effective Hamiltonian which is only 
non- vanishing near the impurity location, at \x\ < £k- In 
a long wavelength effective Hamiltonian this interaction 
appears right at the origin. This interaction is irrelevant 
in the re-normalization group sense; treating it in per- 
turbation theory, its effects become weaker and weaker 
at lower energies and longer lengths. Building on earlier 
work by Wilson^ Nozieresi - developed an elegant de- 
scription of this physics, following the Landau approach 
to Fermi liquid theory, by writing an expression for the 



phase shift which depends on the energy of the electron 
being scattered and also on the local density of electons 
at the origin. In [32j this approach was transcribed into 
an effective Hamiltonian, an approach which is simpler 
to deal with in some situations and more in accord with 
modern renormalization group (RG) techniques. (In a 
similar way, the Landau approach to Fermi liquids can be 
recast in terms of an effective Hamiltonian which can be 
studied by RG methods. See, for example. [33j.) Further- 
more, in [32[ it was shown how this effective Hamiltonian 
could be uniquely determined, up to one overall coupling 
constant with dimensions of inverse energy, 1/Tk, by us- 
ing the one-dimensional (ID) formulation of the Kondo 
model and ID spin-charge separation. Essentially all low 
energy long distance properties of the Kondo model can 
be calculated by perturbation theory in this Fermi liq- 
uid interaction. They are generally proportional to the 
first or second power of the coupling constant, 1/Tk or 
1/T K . It is not an easy matter to determined the precise 
value of Tk for a given microscopic model. However, for 
a weak bare Kondo coupling, v J <C 1 (where v is the den- 
sity of states) it is known that Tk is exponentially small, 
Tk oc D exp[— 1/(JV)], where D is a cut-off or bandwidth 
scale. Furthermore, when the bare coupling is weak, the 
physics is universal at all energy scales <C D. This im- 
plies, in particular, that the ratios of various physical 
quantities calculated in perturbation theory in 1/Tk are 
universal. The most famous such universal ratio, between 
the impurity susceptibility and specific heat, is known as 
the Wilson ratio. However, many more such universal ra- 
tios can be readily calculated. In this section we extend 
extend this approach to the persistent current. 



A. The Boundary Conditions 

Since Fermi liquid theory applies at low energies and 
long distances, it is convenient to linearize the spectrum 
around the Fermi points ±fcp and to introduce left (L) 
and right (R) moving chiral fields: 



(2.1) 



The non-interacting part of the continuum Hamiltonian 
becomes: 



Ho = 



4>iid x ip R - tpXid^i 



(2.2) 



We have set vp = 1. The Kondo interaction becomes: 

H K -> JkS- [4(0) + ^ (0)1^(0) + V>l(0)]. (2.3) 

The strong coupling fixed point of the Kondo problem 
corresponds to the boundary condition: 



c o cxVl(0) + ^(0) = 0. 



(2.4) 



We may think of the sceening electron as sitting at the 
origin, j = and the other electrons (low energy quasi- 
particles) must not enter or leave the origin in order to 
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avoid breaking up the Kondo singlet. Considering the 
model defined on a ring, and ignoring, for the moment, 
the magnetic flux, there is a second boundary condition 
at the strong coupling fixed point: 



At half-filling this becomes: 



cp oc e 



-—ikpL 



^Pr(L) = 0. 



(2.5) 



This strong coupling fixed point corresponds to an open 
chain of L— 1 sites [j = 1, 2, 3, . . . [L— 1)] with no impurity 
spin. Initially we consider half-filling. In this case, for 
such an open chain, kp = tt/2 regardless of the parity 
of L. Since the boundary conditions of Eqs. (|2.41 12. 5p 
are true at all times, while ipR/L is a function of t =p x 
only, it follows from Eq. (|2.4|) that we can regard ipn(x), 
for x > as the analytic continuation of iPl(x) to the 
negative x axis: 

Mx) = -1>l(-x), (x > 0). (2.6) 

Likewise, Eq. (|2.5|) implies that we can also regard ippt(x) 
as the analytic continuation of 

i' R (L-x) =-e- 2lkFL ^ L (L + x). (x>0) (2.7) 

Letting x — > L in Eqs. (|2.7j) . we see that: 

V^(2L) = e 2i ^ L ^(0). (2.8) 

I 



t/> L (2L) = (-1) L M0), 



(2.9) 



periodic for L even and anti-periodic for L odd. We may 
formulate the model in terms of left-movers only with 
these periodic or anti-periodic boundary conditions. 



B. The Fermi Liquid Interaction with Zero Flux 



We now wish to write the Fermi liquid interaction 
in terms of this left-moving formulation of the model. 
This must involve the even channel fermions (even under 
x — ► L — x) only since only the even channel appears 
in the Kondo interaction. Furthermore, it only involves 
this channel near x = 0. Consider for example the even 
channel fermions in the lattice model at a distance of 1 
site from the origin: 



a + c L oc e ifeF Vfl(l) + e- lkF 4> L (0) + e ik ^ L - x ^ R (L - 1) + e^^-V^L - 1). 

I 



(2.10) 



We now use the facts that kp — tt/2, the continuum fields 
iPl(x) and iI>r( x ) are slowly varying on the lattice scale 
and the boundary conditions of Eq. (|2.4l 12. 5p to write 
this in the purely left-moving formulation as: 



ci 



+ c L cx e-"/ 2 V>i(0) + e-W-WfaiL). (2.11) 



We expect only this combination of fields to appear in 
the Fermi liquid interaction, up to higher dimension op- 
erators involving derivatives of the fermion fields. This 
follows from observing, using the boundary conditions of 
Eq. (|2.4[ 12. 5j) and the fact that the continuum fields vary 
slowly, that any non-vanishing even lattice fermion field 
near the origin is proportional to this one. For j even 
and j <C S,k one finds: 

Cj- + c L -j « 0. (2.12) 
While for j odd and j <§C we have: 

+ c L -j oc e- inj ^ L (0) + e -™( L ^)/ 2 ^(L) 
= (-l)W- 1 )/ 2 [e-^M^ + e-^-^ML)] , 

(2.13) 

which is proportional to Eq. (|2.1ip . 



Once we have written the Kondo interaction in terms 
of the even sector fermions only, it is convenient to 
bosonize, introducing separate spin and charge bosons. 
Importantly, the Kondo interaction involves only the spin 
bosons in the even sector. It then follows that the Fermi 
liquid interactions can be written in terms of these bosons 
only. The only possible dimension 2 operator, in the spin 
sector, which respects the SU(2) symmetry is J^(0), the 
square of the spin density operator. Going back to the 
fermion representation: 



J e (0) = ^(0)^ e (0). 



(2.14) 



where, from Eq. (|2.1ip . 



ip e {Q)(x -j= . (2.15) 
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Thus we write the leading irrelevant Fermi liquid inter- 
action, for the case, L even as: 



x [^ L (0)-e-^ L ^ML)]} 2 - (2-16) 

We have written the coupling constant in front of this 
operator as 8tt/3Tk- It follows from standard scaling ar- 
guments that this coupling constant should be of order 
the RG cross-over scale. In general, it should be written 
as 1/Tk times a dimensionless constant of 0(1). All low 
energy propertites of the system can be determined from 
this interaction term, including the impurity susceptibil- 
ity and specific heat, for example. The sign in Eq. (|2.16[) 
is known to be the correct one since it gives the correct 
(positive) sign for these two quantities. The factor of 
87r/3 is purely a matter of a convention, or a precise def- 
inition of what is meant by Tk- Unfortunately, a great 
number of different conventions for Tk are in current 
use. We have chosen here to use the same convention as 
used by Nozieres^ (who also referred to 1/Tk as a) and 
by Glazman and Pustilnik? 4 . In Appendix B we briefly 
review other definitions of Tk in popular use and the 
constant factors relating them to each other. 

Strictly speaking, one other operator of dimension 
2 is permitted, J e (0) 2 , the square of the charge cur- 
rent, rather than the spin current. As mentioned above, 
naively the charge sector decouples from the Kondo in- 
teractions so that this operator should apparently not 
be generated. However, irrelevant operators (at the 
weak coupling fixed point) couple charge and spin sec- 
tors together ultimately permitting this operator to oc- 
cur. However, the coefficient of this operator is expected 
to be 0(1/D), where D is the bandwidth (t in the tight- 
binding model), rather than 0(1/Tk)- In the weak cou- 
pling limit, where Tk <C D, this other interaction can be 
ignored. 

C. The Fermi Liquid Interaction with Finite Flux 

So far we haven't mentioned the flux, $ = ca/e. The 
phases in the hopping terms corresponding to any de- 
sired flux can be inserted anywhere in the ring. They 



can be moved around freely by phase redefinitions of the 
electron fields (gauge transformations). For purposes of 
understanding the strong coupling fixed point it is con- 
venient to imagine that they are inserted far from the 
impurity compared to £k so that the strong coupling 
boundary conditions, discussed above, are unaffected by 
the flux. Of course, this is only possible for L^> ^k but 
it is precisely that limit which we are now considering. 
The Fermi liquid theory only applies in that limit. It 
can be seen that adding phases to hopping terms nec- 
cessarily couples even and odd sectors together, since it 
breaks parity. For this reason, it is more convenient to go 
back to left and right movers and then ultimately to left 
movers only on a ring of length 2L, as discussed above. 
It is also most convenient to put the phase at the origin 
after establishing the strong coupling b.c. This amounts 
to: 

t/j L (L) -> e' la 4, L {L). (2.17) 
Hence the Form liquid interaction becomes: 

x {MO)-e~ l{ * L/2+a) ML)}} 2 - (2-18) 



D. Perturbation Theory 

The remaining calculations are quite straightforward. 
We simply do first order perturbation theory in 1/Tk 
for the groundstate energy, imposing periodic or anti- 
periodic b.c.'s, Eq. (|2.9p . on the purely left-moving 
fermion fields. 

Now we turn to calculating the expectation value of 
H int . From Eq. |[2T8]) . recalling that L = N - 1 at half- 
filling, it is clear that it is advantageous to define a shifted 
flux a in the following manner: 

a = a + Tr(N - l)/2 (N odd) 

a = a + irN/2 (N even). (2.19) 

[This definition can been seen largely as a way of defining 
the coefficient in front of sin(a) in the persistent current 
to be positive.] Using this definition, we first rewrite Hi nt 
as: 



H m t = ^{e 2i& [^(L)a^(0)] 2 -2ie i& ^(L)a^(0)- [^(0)^(0) + ^ (L)<j^(L)] + h.c. + . . .} , (N even) 

H int = ^{-e 2ia [^(L)a^(0)} 2 + 2e i& ^ (L)a^(0) ■ (0)^(0) + ^(L)a^(L)] + h.c. + ...}. (N odd) 

(2.20) 

I 

Here we have dropped the L subscripts and also dropped terms which are independent of a and hence won't con- 
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tribute to the current in first order in 1/Tk- The first 
term can be rewritten: 



H 



(i) 



7r(-l 



-e 2ia ^ Tt (L)Vr(Q)V' ;t (i)V'i(0) + h.c. 



(2.21) 

We now evaluate the various terms using Wick's theo- 
rem. For this we need the equal time propogator for a 
left mover with anti-periodic and periodic boundary con- 
ditions on an interval of length 2L. 



1. N Even at Half-filling 

For the case TV even, with anti-periodic boundary con- 
ditions, this is: 



Now consider the cos a term in the energy. This is 
actually zero for N even but let's go through some steps 
which are useful also for N odd. We use Wick's theorem 
to write this as: 

< ^(L)a<ip(0) • [^(0)^(0) +^(L)aip(L)} > 
= 2 < ^ (£)dty(0) > ■ < ^ (O)oty(O) > 



(2.30) 

For N even it is easy to see that the first term in Eq. 
(I2.30|) vanishes and the last two cancel using Eqs. (|2.23|) - 
(I2~27l) . 



1 OO 



/2)x/L 



-iS% 



4Z sin(7rx/2L) ' 



In particular: 



(2.22) 



(2.23) 



Of course, 

< V tQ (i)#(0) >= - < M0)^ a (L) > . (2.24) 

We also need an expression for < %j}^ a (0)^(0) >. The 
sum in Eq. (I2.22| is ultraviolet divergent so we intro- 
duce a (dimensionless) cut-off, an upper bound on the 
summation variable m: m < _D, giving: 



5<SD 



(2.25) 



We will see that our results for the current do not depend 
on the value of D. It follows from PH symmetry that: 



5<$D 

< ^(0)^(0) >=-§£-• 



(2.26) 



Furthermore 

< M0)^ a (0) >=< ML)^ W (L) >, (2.27) 

as follows from the reflection symmetry of the model 
which takes x — + L — x. 
The cos 2a term is: 



E£> = — —c, — + h.c. = 
T K Ul) 8T K L 2 



With a current: 



j {2) = -edEo/di 



sin 25 (N even). (2.29) 



2. N Odd at Half-filling 

Now consider N odd. The strong coupling groundstate 
has spin -1/2; we choose the state with total S z = 1/2. 
The propogators are now different for spin up or down. 
We find: 

< ^(x)^(0) >= J- J2 ^ mx/L = 77 ( wry 
2L ^-^ AL sm[irx/2L) 

m=o v ' 1 

On the other hand, the J, propagator is the same except 
that the m — term is omitted from the sum, subtracting 
1/2L to the final result. Thus: 



Again we impose a cut-off so: 



Now we have: 



< ^ T (0)^(0) >= 



D — 1 
2L 



(2.32) 



(2.33) 



(2.34) 



The (—1) reflects the breaking of particle-hole (P-H) sym- 
metry. We have one unpaired spin up electron, sitting 
right at the Fermi surface. Furthermore, 



<^t (o)^(o) >= 



D- 1 
2L 



and 



< ^(0)^(0) >= R 



(2.35) 



(2.36) 



cos 2a, (2.28) We may combine these as: 



<^t(0)^(0)> = ^i/^ + -L(tr*)2 
<^ Q (0)V/ t (0)> = ^=^£-^(^2.37) 
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We now substitute this into Eq. (|2.30[) . It is easy to 
see that the terms containing D all cancel, for the same 
reasons that the entire expression vanishes for N even. 
The remaining terms are: 



< ^(L)a^(Q) • [^(0)5^(0) +^(L)ai/j(L)] > 

2 

[2tr(aa z ) ■ tr(aa z ) - 2tr(aa z ■ aa z )\ 



1 

Tl 
i 

Jl 

3 

4L2- 



[8 + 4] 



Thus the cos a term in the energy, for N odd, is: 



E, 



(i) _ 



The term in the current is 
.(i) _ en 



.2e i& —+cc - 77 
~ 4I/ 2 ~~ 2TkL 2 



cos a. 



J 



2T K L* 



since (N odd). 



(2.38) 



(2.39) 



(2.40) 



The sin(2a) in the current for odd N is given again by 
Eq. (|2.29|) . with the opposite sign: 



;(2) _ 



4Tffi 2 



sin 2a (N odd). 



(2.41) 



Note that the minus sign which follows from spin up and 
down Green's functions having opposite sign for N odd, 
Eq. (|2.32p . replaces the minus sign due to the factors of 
i in the Green's function for TV even, Eq. (|2.23p , leaving 
only the minus sign from Eq. (|2.20|) . 

We have so far set vp = 1. We summarize our results 
on the persistent current, for <C L below, reinserting a 
factor of vp by dimensional analysis and replacing vp /Tk 
by S,k ■ 



Je 
jo 



ev F £ K ir 
L L 4 

evF £k 
L L 



sin 2a 



sin a 



■ sin 2a 



(2.42) 



It was argued earliei— , that both j e and j can be written 
as evp/L times universal scaling functions of £k/L and 
a. Eq. (|2.42|) indeed has that form. 



3. Off Half-filling 

Now consider the case away from 1 /2-filling, where P-H 
symmetry is broken. Here we refer to particle-hole sym- 
metry breaking by an amount of 0(1) not an amount of 
0(1/ L) as occurs, even at 1/2-filling for N even. Note 



that we continue to choose the reduced band to be sym- 
metric around kp: kp — A < k < kp + A. However, 
the process of integrating out wave- vectors further away 
from kp is not P-H symmetric, when kp does not have 
the P-H symmetric value, so we expect to generate PH 
symmetry breaking terms in Hint. The most important 
such term, which can lead to a current of 0(1/ L), is: 

H 2 = -e" (A, - 1 )/ 2 Ae M ^ t (i)V'(0) + h.c. (2.43) 

where A is a real coupling constant. Here the phase in this 
expression is determined by parity, i.e. by the fact that it 
arises from a term —2\if>\^) e using Eq. (|2.11[) . This term 
only occurs when P-H symmetry is broken. We expect 
A oc (Jx/t) 2 where t is the original bandwidth. 

More generally, we might start with a microscopic 
model which includes second neighbor tunnelling past the 
quantum dot: 



SH = -t 2 (4_ 1 ci +h.c). 



(2.44) 



Such a term breaks P-H symmetry even at half-filling and 
leads to a term in the low energy effective Hamiltonian 
of the form of H 2 in Eq. (|2.43|) with \ott 2 . 

If the system is very close to half-filling, n — 1 = 8n <C 
and we begin with a pure Kondo model with no 
potential scattering or direct tunnelling, then we can cal- 
culate the value of A, in H 2 directly from the Fermi liquid 
interaction Hi nt of Eq. (|2.20p . To do this, it is convenient 
to integrate out wave-vectors in the RG transformation 
symmetrically around k — -k/2 even thought this is no 
longer fci?, which has the value: 



kp = (tt/2)(1 + 5n). 



(2.45) 



In this way, A is not generated during the RG transforma- 
tion, and we obtain only the interaction Hint- However, 
the low energy effective theory inherits an asymmetric 
cut-off: 

kp-A - (tt/2)Sji <k<k F - (7r/2)5n + A. (2.46) 

We may now generate H 2 from Hi nt by simply normal 
ordering Hi nt . This follows since Eqs. (|2.25l I2.26P are 
now modified to: 





' D 


Sn 


°H 


2L 


+ T 


SOL 


' D 


Sn 


6 P 


2L 


~ T 



(2.47) 



We now write H int as a normal ordered part plus a cor- 
rection of the form of H 2 . This follows using: 
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^(L)aip(0) • [^(0)^(0) +^(L)a^(L)] =: ^(L)aip(0) • [^(0)^(0) + tf(L)ffij){L)] 
[^(L)M0) < 0|#(0)V 7t (0)|0 > -^ 7t (i)#(0) < 0\^(L)ML)] (<?f ■ d*). 

I 



(2- 



Using Eq. (|2.47[) we see that the two terms in the second yielding instead: 
line of Eq. (|2.48p don't cancel, away from half-filling, 



of 

ipi(L)aip(0) • [^(0)^(0) + ^(L)aifj(L)] = : ^(L)ai/j(0) • [^(0)^(0) + tp j {L)aip(L)] : (L)i/j(0). (2-49) 



r 



Inserting this expression into Eq. 
quadratic term of the form of H2 with: 



A = 



nSn 
27V 



(|2~20f gives a 



(2.50) 



for N even. 

We may now evaluate the additional terms in the cur- 
rent arising from H2, using perturbation theory in Hint- 
We only discuss first order perturbation theory. We eval- 
uate < H 2 > using Eq. (I2.23p . This gives an additional 
current: 



Sj = (eX/L) sin a, 



(2.51) 



for N even. For N odd, there is no extra contribution to 
first order in A since < ip> a (L)ipp(0) >cx (cr z )p- Hence, 
to first order, we expect that in the limit of small Sn the 
current for N odd, j a , remain unchanged with respect to 
it's value at half-filling. 

Note however, that we have only evaluated the contri- 
butions of Hi n t to the current in first order perturbation 
theory in H2. Thus we haven't ruled out the possibil- 
ity of a term of 0(1/ L) in the current for N odd, with 
broken PH symmetry. 

Adding the universal result for the PH symmetric case 
to this correction from PH symmetry breaking (and rein- 
serting the factor of vf previously set to one) gives a 
persistent current, for N even: 



3t 



evp 



4L 



sin(2a) + Asin(a) 



(2.52) 



The non-universal term arising from PH symmetry 
breaking, proportional to A will always dominate at suf- 
ficiently large L. However, provided that the dimension- 
less coupling, A <C 1, as we expect for small bare Kondo 
coupling and not too strong direct tunnelling across the 
quantum dot, the first, universal, term will dominate over 
the range of lengths: 



6c < L < 6f/A. 



(2.53) 



The presence of a term in j e proportional to sin(<5) in 
the absence of PH symmetry was shown in Ref. |30j using 
strong coupling perturbation theory. The above deriva- 
tion provides further independent analytical evidence for 
such a term. 



III. NUMERICAL RESULTS 

The Fermi liquid theory results for the persistent cur- 
rent developed in the previous section are parameterized 
in terms of Tk and should be valid in the regime T <C Tk, 
1 <C £,k ^ L. For a useful numerical test of these results 
it is therefore necessary to study ratios of the Fourier 
components of the persistent current that become inde- 
pendent of Tjc and can be tested at T — using exact 
diagonalization (ED) methods. Furthermore, in order 
for 1 <C £k to hold, Jk should not be too large. A more 
severe constraint is that the system sizes should satisfy 
L 3> £r- • With an exponentially diverging ^k at small 
Jk we see that we quickly leave the regime of validity 
of Fermi liquid theory as Jk for the rather modest 
systems sizes we can treat using ED. Clearly, the inter- 
esting regime is then intermediate values for Jk- Our 
approach is then to calculate the different Fourier com- 
ponents, a n , of the persistent current, j, as a function of 
the flux, a. Our convention for the Fourier components 
are j(a) — J2 n a « sin(na). Note that, a very precise de- 
termination of the a dependence is necessary in order to 
determine the Fourier components reliably and we typ- 
ically use 200 values for a. The numerical results can 
then be compared to the predictions of the Fermi liquid 
theory as well as results previously obtained using strong 
coupling perturbation. 30 



A. At Half-filling 

We start by discussing our numerical ED results ob- 
tained at half-filling with systems sizes of L = 12 — 15. 
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1/J K 1/J K 



FIG. 1: The ratio, a^/al, of the two first Fourier coefficients of 
the current for the SCQD at f /2-filling, for TV odd. o indicates 
ED results with L — 12 and N = 13 while • indicates ED 
results with L — 14 and N = 15. The dashed line (barely 
visible) is the strong coupling perturbative result, Eq. (|3.2|1 . 
for this ratio with L — 12, the solid line is L — > oo the strong 
coupling perturbative result, Eq. (|3.3|l . the dashed dotted line 
the fermi liquid result, Eq. (|3.4[1 . and the dotted line the ideal 
ring result, Eq. (|3.5|l . 



FIG. 2: The ratio, a^/a^, of the second Fourier coefficients 
of the current for the SCQD at 1/2-filling, for TV odd and 
TV even, o indicates ED results with L = 11, TV — 12 and 
L — 12, TV = 13 while • indicates ED results with L — 13, 
TV = 14 and L = 14, N = 15. The dashed line is the 
strong coupling perturbative result, Eq. (|3.2[) . for this ratio 
with L cvcn = 11, i odd = 12, the dashed dotted line the fermi 
liquid result, Eq. (|3.4p . and the dotted line the ideal ring re- 
sult, Eq. ([33J. 



Due to the lack of symmetry in these systems it is dif- 
ficult to reach larger sizes and the largest system (L = 
14, TV — 15) required a diagonalization of a matrix of size 
20,796,633 for each value of a after symmetry reduc- 
tions. 

In Ref l3fj strong coupling perturbative results were 
given for the persistent current at half-filling past the 
SCQD. It was shown that, for odd or even TV, 



joL/e 



jeL/e 



32 . 7T . . 37T . 

— [ tan (— ) + tan (7rr)J sma 



9 J 



V 2L ; 



'277 



K 

12 8 r„ . - . ,„~v, 

'iJ^L 1 sma - sm ( 2a )] 
32 



3J K L 



[1 + l/cos(7r/L)] 2 sin25. (3.1) 



As elsewhere, L does not include the impurity site while 
TV does include the impurity electron. The definition of a 
is also the same: For TV odd a = a + ir(N — l)/2 and for 
TV even a — a + irN/2. If we concentrate on the first two 
Fourier coefficients of the current we see from the above 
strong coupling results of Eq. ()3.1|) that: 



-12 



J K L[tan(7r/2L) + tan(3?r/2L)] 

-4 L e 
(l + l/cos(7r/X e )) 2 Z°' 



24 



(3.2) 



Here the superscripts o and e refer to N odd or even, 
respectively. In the limit L — > oo we see that the strong 
coupling expansion gives: 



-6 



ttJk + 12 



^ -> -1. 



(3.3) 



The Fermi liquid theory developed in the preceeding sec- 
tion, valid for 1 < ^- < L, predicts, using Eqs. (|2.29|) 
and (|2~10"|) : 



-1. 



(3.4) 



When Jk is exactly zero we obtain the behavior of an 
ideal ring, in this case the current has a charateristic saw- 
tooth like for m 25 i 35 where the same Fourier coefficients 
can be trivially calculated. One finds at J# = : 



o|_l o|_ 



(3.5) 



Note that this result is dramatically different from the 
result of Eq. I|3.4p . valid for 1 <S £ K < L. The ratio of 
the Fourier coefficients appear discontinuous in the limit 
Jk -» 0. 

We begin by discussing our results for a^/aj which are 
shown in Fig.Q]for system sizes of L — 12 and L — 14, for 
a range of different Jk- The theoretical strong coupling 
result of Eq. (|3.2p for L — 12 is shown as the dashed 
line and differs only slightly from the L — > oo result of 
Eq. (|3.3[) shown as the solid line. The numerical results 
agrees with the strong coupling results once Jk > 10. 
In the opposite limit, Jk — ► 0, we see that the numerical 
results quickly approach the ideal ring result of Eq. (|3.5|) . 
shown as the dotted line. This is very reasonable, since 
the very limited system sizes of L = 12, 14 certainly no 
longer satisfies the equality L 3> for Fermi liquid 
theory to be valid once £k diverges as 0. The 

interesting region is therefore the intermediate coupling 
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region, 0.1 < Jk < 10. A very rapid crossover from the 
ideal ring result of 1/2, at small Jk, to a negative value 
is observed, consistent with an exponentially diverging 
Furthermore, a^/a" develops a pronounced plateau 
("dip") around Jk ~ 4, reaching negative values, even 
for these rather modest system sizes. As we increase the 
system size from L = 12 to L = 14 this "dip" becomes 
slightly more pronounced. This implies that the cross- 
over between the ideal ring result for this ratio of 1/2, 
at small Jk, to the Jk — > oo result of 0, can not be 
monotonic even in the thermodynamic limit. We take 
this to be indicative of the validity of the Fermi Liq- 
uid theory. We expect that numerical results for a%/a° 
in the thermodynamic limit would roughly follow the 
strong coupling result given by Eq. (|3.3|) (the solid line 
in Fig. [T]) out to Jk ~ 0.1 and then attain the value -1/2 
rather quickly, jumping discontinuously for Jk = 0, in 
accordance with the Fermi Liquid theory. In the subse- 
quent section, where we discuss our results away from 
half-filling, we therefore exclusively focus on the value of 
Jk = 4 where the "dip" occurs, since this would appear 
to be the most promising value for Jk for observing FLT 
behavior with the available system sizes. 

Results for a^/a% are shown in Fig. [21 In this case 
the behavior as a function of Jk is monotonic and the 
evidence for a regime described by FLT perhaps less ob- 
vious. The dotted, dashed and dashed-dotted lines are 
the ideal ring results of Eq. (|3.5|) . the strong coupling re- 
sult of Eq. ([HI]) and the FLT result of Eq. ((S3). Again 
an extremely rapid crossover is seen at intermediate val- 
ues of Jk- Furthermore, as the system size is increased 
from L = 12 to L = 14 this crossover moves to smaller 
values of Jk approaching the predictions of Fermi Liquid 
theory. 



B. Off Half-filling 

Next, we turn to a discussion of our results for off half- 
filling. With our definitions of N and L we define, as 
elsewhere, Sn — (N — 1)/L — 1. A typical Hilbert space 
size is 64,414,350 for [L = 15, N = 14) after symmetry 
reductions. In Ref 30 strong coupling perturbative results 
for the current past the SCQD away from half-filling were 
also given. At large Jk it was found that for N odd and 
N even: 



j L/e i 



32 
9J$ 



sin(<S) 



Tr(iV-l) 7T 3ttOV-1) 3tt 

sm tan sm — tan — 

2L 2L 2L 2L 



JeL/e 



32 
9jJ 



sin(a) 



cos(ir{N - l)/(2£)) cos(3tt(7V - l)/(2£)) 



cos(tt/(2L)) 



cos(3tt/(2Z)) 



CO 



0.20 



0.15 



r0.10 



0.05 



0.00 



' — I — ' — ' — 

• • ED, L=13 
▲-▲ ED, L=14 

m-m ed, L=15 



-0.8 



-0.6 



Sn 



-0.4 



-0.2 



FIG. 3: The first Fourier coefficient scaled by L, afL, of the 
current for the SCQD away from 1/2-filling, for N even as a 
function of Sn = (N — 1)/L — 1. In all cases Jk = 4.0. o in- 
dicates ED results with L = 13 and N = 14, 12, 10, 8, 6, 4, A 
indicates ED results with L = 14 and N = 14, 12, 10, 8, 6, 4, 
finally • indicates ED results with L = 15 and N = 
16, 14, 12, 10, 8, 6. The lines are guides to the eye. 



+ o, 

Eq. 



joL/e 



-(Sn) 2 sin(5)/( J 2 K L). 



This results predicts that as L — ► oo, Sn - 
is changed from the value at half-filling, 
only by a term proportional to 

If we only work to linear order in Sn this correction 
can be neglected. However in this limit, j e L/e ~ 
sin(a)647T(5n/9 J K , that is, a sin(ci) term proportional to 
Sn. At half-filling this term is clearly for j e , and the 
sin(2a) term in Eq. (|3.ip dominates. We note that this 
limit of the strong coupling results agrees with the pre- 
viously developed Fermi liquid theory that, in the limit 
Sn —>■ 0, predicted no change in j a but the generation 
of a sin(a) term proportional to Sn/ L for j e , Eq. (|2.50[> . 
(I2.51[) . In the following, we therefore focus exclusively 
on the first Fourier coefficient, a±, determining its de- 
pendence on Sn. 

Our results for a\L for ./V even are shown in Fig. [3] 
for system sizes L — 13, 14, 15 for a range of Sn. All 
results are for an intermediate coupling of Jk = 4, the 
most promising coupling to show clear indications of FLT 
behavior for our limited system sizes. Two conclusions 
are immediately evident; since the results fall on a sin- 
gle curve this term in the current is indeed proportional 
to 1/L as predicted by theory. Secondly, for small Sn, 
a\L increases approximately linearly with Sn, again con- 
sistent with the theoretical predictions. At larger Sn the 
results strongly deviate from linear behavior as one would 
expect. 

Finally, we show results for a°L 2 for N odd in Fig. 0] 
for system sizes of L — 12, 13, 14, for a range of Sn. The 
results are again for an intermediate coupling of Jk = 4. 
As was the case for a\L, the results for different system 
sizes again follow a single curve, validating the scaling 
al ~ 1/L 2 predicted by theory. A term proportional to 
(3.6) 1/L can definitely be excluded for this Fourier coefficient. 
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• • ED, L=12 

a-a ED, L=13 
m-m ED, L=14 
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-0.2 



FIG. 4: The first Fourier coefficient scaied by L 2 , a°L 2 , of the 
current for the SCQD away from 1/2-filling, for N odd as a 
function of 8n = (N — 1)/L — 1. In all cases Jk = 4.0. o indi- 
cates ED results with L — 12 and N = 13, 11, 9, 7, 5, 3, A indi- 
cates ED results with L = 13 and N = 13, 11, 9, 7, 5, 3, finally 
• indicates ED results with L = 14 and N = 15, 13, 11, 9, 7, 5. 
The lines are guides to the eye. 



first term is referred to as the inelastic part, corresond- 
ing to an electron-electron interaction at the origin. (For 
a derivation of the Fermi liquid interaction in this form, 
see [32[ , but note that the fermion fields are defined with 
an unconventional normalization there so that they are 
larger by a factor of V2n.) Now the phase a, is regarded 
as a time-dependent vector potential, related to the po- 
tential energy drop across the quantum dot by: 



o+ 



da/dt = - / dxE{x) = AV. 



(4.4) 



(In this section we set the electron charge, e = 1, reinstat- 
ing it at the end.) The conductance through a quantum 
dot, treating the leads as ballistic, is insensitive to where 
the electric field is applied so we have chosen, for conve- 
nience, to apply it right at the junction, over a region of 
width of order 

We write the current operator as the rate of change of 
one-half the difference of the total number of electrons 
on the right hand side (x > 0) of the system minus the 
total number on the left hand side: 



We also see that, for Sn — » 0, our results are consistent 
with deviations from the result at half-filling being of 
higher than linear order in Sn, although a definite con- 
clusion is hard to obtain due to the limited system sizes 
available. 



IV. FINITE TEMPERATURE CONDUCTANCE 

We now calculate the conductance for a quantum dot, 
side-coupled to infinitely long leads, at a low finite tem- 
perature, T « Tk- We do this using our Fermi liquid 
Hamiltonian, of Eq. (|2.18[) . We choose L and L/2 even, 
for convenience but we are now taking the L — > oo limit 
so this choice is immaterial. Working with left-movers 
only, and taking into account that the point x = L corre- 
sponds to x = 0~, x = to + , we write the Fermi liquid 
interaction as: 



Hint 

Here 



3T K yv 2 Vi 3T K 



rp^ip ■ -(-const. (4.1) 



i>= -j= [^ L (0+)-V^(0-)e- M ] , (4.2) 



the : : denotes normal ordering and we define: 



f{x)£g(x) = f 



dg_ _df_ 

dx dx 



a- 



(4.3) 



The second, derivative term in Eq. (|4.1[) . arises from a 
point splitting procedure. The second term is referred 
to as the elastic part of the interaction, corresponding to 
a single electron impuirity scattering process, while the 



I = (l/2)d(N+ - N-)/dt = (l/2)i[H,N+- N-], (4.5) 
where 

/>OC 

N+= [rt>tix)M*)+^x)M*)]> ( 4 - 6 ) 



and similarly for N- . As discussed above, we can use the 
perfectly reflecting strong-coupling b.c. of Eqs. (|2.4p . 
(12. 5|) ) to rewrite the left and right movers at x > in 
terms of left-movers only on the entire real line. Thus we 
define: 



4>l+(x) = iJjl(x), (x > 0) 
= ii R {-x) (x < 0). 



(4.7) 



Here the subscript + indicates that this field derives from 
the x > region. Similarly we define another left-moving 
field on the entire real line, ^Pl~{x) which derives from 
the region x < 0: 



i>L-(x) 



In this notation: 



N+ = 



i/> L (x), (x<0) 
^r(-x) (x > 0). 



ipl ± (x)ip L± (x) 



(4.8) 



(4.9) 



and the operator appearing in the Fermi liquid interac- 
tion, defined in Eq. (|4.2[) becomes: 



1 



V2 



[^+(0)-Vi-(0)e 



(4.10) 



We wish to calculate < I > in linear response to a 
time dependent vector potential a(t). At this point, there 
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are two ways of proceeding. 
I = (l/2)d(AN)/dt where 

AN = N, 



We may use the expression 



(4.11) 



and ultimately relate the conductance to the Green's 
function of AN, or we may calculate explicitly the com- 
mutator [H, AN] and relate the conductance to the 
Green's function of that operator. In both approaches 
the calculation must be carried out to second order in 
the Fermi liquid interaction and the amount of work is 
roughly the same either way. However, it is actually 
much more convenient to work with dAN / dt because the 
needed Green's function can be related to the T-matrix, 
which was calculated previously 10 ' 32 This approach was 
used in [34[ for instance, to calculate the conductance 
through an embedded quantum dot. The calculation is 
very similar in the side-coupled case and we may obtain 
the answer by only a slight modification of their result. 
This approach is used in the rest of this section. It is 
also of some interest to calculate the conductance using 
the commutator method and this is done in Appendix A, 
yielding precisely the same result. 

In the dAN /dt method, we use the fact that, to 0(a), 
H = H(0) + a(t)I, to obtain the conductance: 

C = lim - / dte luJt < [I(t),I(0)} > . (4.12) 

w^O LO J Q 

Integrating by parts: 

/>oo 

C = (1/4) lim u / dwe iut < [AN(t), AN(Q)} >, 
Jo 

(4-13) 

where the retarded Green's function is evaluated at a = 
0. It is now convenient to go to the even-odd basis: 

V'e/o = J=Wi+T^-). (4.14) 

Since the interactions only involve ip e , C factorizes into a 
free Green's function for ip a multiplied by the non-trivial 
Green's function for ip e . This latter Green's function can 
be expressed in terms of the T-matrix, giving a formula 
for the conductance^ 



C= — 



E 



de(-df/de)[-KuImT s (e)}, (4.15) 



where s = ±1 labels the fermion spin, /(e) is the Fermi 
distribution function at temperature T and v is the den- 
sity of states. As shown in [34| following [13, for the 
embedded quantum dot: 



-wvT s (e) 



where 



i [exp[2iS s (e)] 
ii 



l]+cxp[2z(5 s (e)][-7n/n„(e)], 
(4.16) 



6(e) 
-irvf in (e) 



sir /2 + 6(e) 
u>/T K 

, (e 2 + 7T 2 T 2 ) 
* 2T 2 

K 



(4.17) 



6(e) and T,- n (e) are the phase shift and inelastic part of 
the T-matrix calculated in perturbation theory in the 
Fermi liquid interactions. The extra stt/2 term in 5 S 
arises from the ±7r/2 phase shift characterizing the strong 
coupling fixed point. This ±7r/2 phase shift reflects the 
perfect transmission (C = 2e 2 /h) at T = for the em- 
bedded quantum dot. On the other hand, for the side- 
coupled quantum dot, C = at the zero temperature 
fixed point, implying the absence of this extra phase shift. 
Thus, for the side-coupled quantum dot we have: 



-irvT s (e) 



1 

27 



1 



+exp[2^ s (e)][-7ri/Tn(e)]. 



exp [2iS s (e)} 

. (4 ' 18) 

To order u> 2 , T 2 , the difference between T- matrices for 
the embedded and side-coupled quantum dots is just a 
change in sign in the first and third terms. This implies: 

-TwImT(e) = 1 - ^ J~ 2 T ^ (embedded) 



2T 2 



-7TfImT(e) 



(3e 2 + n 2 T 2 ) 



2T 2 

K 



(side-coupled) (4.19) 



Inserting these expressions into Eq. (14. 1 5[) for the con- 
ductance, and restoring h and e which were previously 
set to one, gives: 



G 



h 


1 - 


2e 2 tt 2 T 2 


h 


T 2 

K 



7T 2 T 2 



T 2 



(embedded) 



(side-coupled) 



(4.20) 



We may now write universal Wilson-type ratios be- 
tween the low temperature conductance, C, through a 
side-coupled quantum dot and the zero temperature per- 
sistent current, j e , j a , through the same quantum dot 
inserted into a ring of size L ^> £k '■ 



L 2 je/ev 2 F 
yJhC/(2e 2 T 2 ) 

L 2 jo/ev 2 F 
^hC/(2e 2 T 2 ) 



— sin 2a 

4 

1 . - 1 . „- 

— ana sin 2a. 

2 4 



(4.21) 



When PH symmetry is broken the conductance also 
gets a contribution of second order in A the coupling 
constant in the term, H2 of the effective Hamiltonian, 
given in Eq. (|2.43p . The simplest way of evaluating this 
contribution is to observe that, since the Hamiltonian 
is non-interacting, ignoring the other term Hi nt in Eq. 
(|2.18p . we can evaluate the conductance at zero temper- 
ature using the Landauer formula, C — (2e 2 /h)T r where 
T r is the transmission probability through the quantum 
dot. For small A, T r = A 2 . 

Including both terms, the conductance takes the form 
at low temperatures: 



C 



tt 2 T 2 
T 2 



+ A 2 



2f_ 
h ■ 



(4.22) 
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At sufficiently low T the second, non-universal term, aris- 
ing from PH symmetry breaking always dominates. How- 
ever provided that the PH symmetry breaking (and, in 
particular, the direct tunnelling across the quantum dot) 
is small, the first, universal, term dominates for: 

XT K « T « T K . (4.23) 

Note that this situation is very analogous to that for the 
persistent current, with T replaced by Vp/L. 

V. CONCLUSION 

We have developed a Fermi liquid theory of the persis- 
tent current in a side coupled quantum dot. This theory 
should correctly describe the limit a <C £,k *C L where 
a is a microscopic length scale. Numerical results are 
largely in agreement with the existence of a regime cor- 
rectly described by this Fermi liquid picture. However, 
probably due to the limited system sizes available and 
the requirement that L 3> S,k, the numerical evidence 
cannot be described as strong. Our Fermi liquid theory 
confirms the existence of a term in the persistent current 
proportional to sm(a)Sn/ J%L for N even, absent at half- 
filling. We have also calculated the conductance through 



a side-coupled quantum dot at low T <C TV an d calcu- 
lated a universal "Wilson" ratio relating the conductance 
to the persistent current. 



APPENDIX A: CONDUCTANCE BY 
COMMUTATOR METHOD 

In this appendix we repeat the calculation of the con- 
ductance, writing the current as I = i[H, AN)/2. This 
provides a check on the previous Fermi liquid calculations 
involving the T-matrix in [To.32.34,36| and also provides 
an instructive example of how universal information can 
be extracted from a cut-off dependent result. 

Since [N+ + N-,H] = 0, we may equivalently use 
/ = i[H,N + ]. Note that N + commutes with the non- 
interacting part of the Hamiltonian which is simply that 
of free fermions with perfectly reflecting b.c.'s at the ori- 
gin, since this Hamiltonian does not transmit any elec- 
trons between left and right sides of the system. The 
commutator of N+ with Hi nt is readily calculated. Since 
all fields in the commutator are left-movers sitting at 
x = 0, we drop the L subscript and the (0) argument for 
simplicity. The result is: 



I = i[H int ,N+] = -Jji{[ e — jpl^_-h.c.],[^+ 



3T K 
1 

2Tk 



V>1 



dx 



ipl^e- ia + h.c 



(Al) 



We note that this is, equivalently, / = dHi nt /da. We 
expand the current operator up to first order in the vector 
potential, a: 



I = Iq + ali + 0(a 2 ). 



(A2) 



Here: 

Io = 

where 



4iri 

3T^ V+ 2 
1 



(7 



2T K 



J+ 2 
h.c. 



nn lei i 



(A3) 



(A4) 



and we have defined the inelastic and elastic terms in 
the current operator which are, respectively, quartic and 
quadratic in fermion operators. The precise form of I\ 
will not be needed. The conductivity is obtained from 
< I >, calculated to first order in a. Thus we obtain: 



I(uj) = [Gr(u)+ < h >}a(uj), 



(A5) 



where Gr is the retarded Green's function: 

G R (w) = -i \ e wt < [I (t), J o (0)] > T ■ (A6) 
Jo 

It can be seen that I — when a(t) is a constant, inde- 
pendent of t. A non-zero / in this case would correspond 
to a persistent current, but this must vanish for infinite 
L since a constant phase a in Hi nt can be eliminated by 
a gauge transformation in that limit. Thus we conclude 
that < h >= -Gfl(0), so that: 

I(u) - [Gr(w) - G fl (0)]aH. (A7) 

From Eq. (]4.4[) . we see that: 

J(w) = C(lu)/\V{uj), (A8) 

where the conductance, C(w), is given by: 



cv) = -[GflH - g r ( 



(A9) 
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where w„ = 2irnT and (3 = 1/T. (We set Boltzmann's 
constant equal to one.) The dc conductance is: 



C= lim -\G r (lj) - G R (0)}. 



(A10) 



We calculate the retarded Green's function of Iq by an- 
alytic continuation from the imaginary time, Matsubara 
Green's function: 



Q{ioj n ) 



P 



dre^ T < / (r)/o(0) > . (All) 



The free fermion Matsubara Green's function, at late 
times, using the same normalization as in Eq. (12.22p . 
that was used in the calculation of the persistent current 
is: 



<v! q (t>mo) >- 



626. 



7 U V 



2/3sin(7rr//3) 



(A12) 



where i, j = ± labels right and left sides of the junction, 
x = ± and /? is the inverse temperature. Using Wick's 
theorem, the identity: 

Y^[tra a a b tro a o h - tro a o h o a o h \ = 18, (A13) 

a, 6 

and collecting the various terms, we obtain the long-time 
behavior of the Matsubara Green's function: 



< Iin(r)Iin(0) >- 



16tt 2 



1 



T 2 K (2/3)4 sin 4 (7rr//3) 



(A14) 



To evaluate < 7 e z(r)/ e i(0) >, we use the derivatives of 
the free fermion Green's function: 



1 



dx \sm[ir(r + ix)/P] 
d 2 1 



i7rcos[7r(r + ix)//3] 
' (3 sin 2 [it {t + ix)//3}' 
it 2 1 + cos 2 [tt(t + ix)/(3] 



dx 2 sin[vr(r + ix)//3] (3 2 sin 3 [7r(r + ix)/0\ 

(A15) 

The terms proportional to cos 2 (ttt//3) cancel leaving: 

8tt 2 1 



< I el {r)I el {0) > 



Tf (2/3)4 sin 4 (^r//3) 

(1/2) < 7 in (r)7 in (0) > .(A16) 



We now consider the Fourier transform of G(r) at low 
frequencies and temperatures. The analytic continuation 
of Q{iui n ) [bj n = 2-Kn/fi) to real frequencies gives the 
retarded Green's function: 

g(iu n -> u + iS) = G r (lj), (A17) 

where 6 is a postive infinitesimal. We will find that 

G(iu n )-> A + B\w n \ + 0(u>l), (A18) 





Mb 



FIG. 5: Contour used to evaluate the integral in Eq. (|A21[l , 



for lo„ << D. The analytic continuation gives: 



dr 



2m 



(-iu) 



dr e (lS -^ 7 



-i(S—iaj)r 



2ni 



(A19) 



G(u) n ), which is an even function of lu h , in our assumed 
particle-hole symmetric model, contains non-universal ul- 
traviolet cut-off dependent terms. Introducing a cut-off, 
To with dimensions of time, and of order 1/D, we find, 
at low u) n : 



1 

T 2 
L K 



+ b 



roP 2 



fl 2 



■ d— + e\ui n 
to 



(A20) 



where a, 6, c, d, e are dimensionless numbers. While the 
terms of 0(w° ) and (lu 2 ) are non-universal and cut-off de- 
pendent, we expect that the terms oc \uj n \ and |<x>„| 3 are 
not. This is related to the fact that only these terms are 
singular functions of u> n . This singularity arises from the 
universal long-time behavior of G(r). Note that the dc 
conductance is completely determined by the universal 
c-term oc \u> n \ and so is a universal quantity in the sense 
that it is independent of the details of the ultraviolet 
cut-off. To verify these assertions, we consider a partic- 
ularly simple ultraviolet cut-off. The t integral defining 
the Fourier transform of the current Green's function is 
restricted to tq < t < [3 — tq Thus we begin with: 



Q{iuj n ) 



-2Att 2 



71(2/3)4 



P-TO 



dre 1 ' 



sin (ttt/P) 



(A21) 



This integral can be straightforwardly developed in an 
expansion in To//3, which has the form of Eq. (|A20[) . This 
can be conveniently done by deforming the r-integral into 
the complex plane. The original integral is equal to the 
sum of an integral along the straight line from r = Itq 
to r = Itq + (3 plus the integral on two quarter-circle 
contours from t — tq to t = itq and from (3 + itq to 
r = [3 — t , which we label C\ and Ci respectively. [See 
fig. ©.] The integral along the straight line can be 
written: 



K 



P d re iu n (r+iT ) 

o sin 4 [7r(r + ir )//3] 



(A22) 
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■iro)/P 



We may now Taylor expand the denominator: 

- oc 
— ifi„47ri(T+iT )//3 2ttm 

sin [tt(t + iT O )/0\ ^ 

(A23) 

Note that the sum converges due to the factors e~ 27rr °/' 3 . 
Here the the binomial coefhcients: 

a m = (m + l)(m + 2)(m + 3)/6. (A24) 

We now integrate over /3 term by term. The result is: 

l dre ,[^ + (2 +m) ( 2 ^)]r =/3(5n _ (2+m) . (A25) 

(Recall that w n = 2im/ (3.) Thus we obtain: 



The integration contour for z is a circle of radius 
e -2jrT //3 < i. This contour encloses a pole of order | n | — 1 
when n < —1. 

We also must perform the integration over the two 
quarter-circle contours: 



dz- 



~{Jc t sin (ttz/P) 



(A27) 



On C\ we write z = r^e 19 where 9 goes from — * ir/2 
along the contour. On C% we write z = j3 + Toe 16 where 
now goes from n/2 — > n. Since the integrand is invari- 
ant under a translation: z — > z + /3, we may combine 
these two terms to write: 



K = - 



3tt 



2ir 



{-u n ). (A26) 



K 12 = ir 



^ e m exp[w w Toe' p J 
o sin 4 (7rr e l V/3) 



(A28) 



Alternatively, the integral of Eq. (|A22[) can be evaluated 
by introducing the complex variable z = e l2 ^( T + lT o)/0 _ 



Next we expand exp[io; rl roe l6 '] in powers of w„to and ex- 
pand sin _4 (7rToe l£, //3) in powers of tq//3, giving: 



K 12 = ir Q I / < i 



i8 -Ai8 



1 + iuj n r e w - 



•2i8 



iuJ n T c 3i8 



r 



'0 „4i6> 



6/3 2 



120/3 4 



(A29) 



We Taylor expand, collect terms, and integrate term by 
term using: 

/>7T 

dde me = 7T (m = 0) 

= (m even, m > 2) 

= 2i/m, (m odd). (A30) 

Note that non-zero terms from the expansion must be 
proportional to e or else e 21 to give a non-zero con- 
tribution. The first type of terms are odd functions of 
u> n , proportional to u> n and w 3 . All other terms are even 
functions of u> n . Keeping only the terms which are non- 
vanishing as To — * gives: 



-2u n f3 2 ujIP 4 



3 TV 



2[3 4 
6tt 3 ' 37r 4 r n 3 



7T 4 T 



4/? 2 
3ir 2 T 



The omitted terms are all even functions of uj n and vanish 
as r — > 0. 

Note that when we add K + K\i we regain an even 



function of u> n , as we must: 



K + K 12 = 



-2/3 2 

2/3 4 
3^ 4 r n 3 



7r 4 r 



+ 



fan 

2ir 

4/3 2 

37T 2 To 



(A32) 



The first two terms are singular functions of co n at w n = 0. 
All remaining terms are non-singular, even powers of w n . 
The singular terms are cut-off independent, unlike the 
non-singular ones. 

Thus we may write the imaginary frequency Green's 
function for the current operator as: 



5(0) 



7rr 2 K| 

T 2 



0{wl 



(A33) 



The analytic continuation to real frequency is straight- 
forward, using Eq. (j A19|) : 



G R (u) Gfl(0) 



iirT 2 uo 

r 2 



0(iu 2 



(A34) 



Thus, we obtain the dc conductance from Eq. (|A10[) 
ttT 2 2e 2 n 2 T 2 



C 



T 2 



h T 2 K 



(A35) 



the same result obtained from the derivative method, Eq. 
(|4.20[) . (A factor of e 2 /h, previously set equal to one was 
inserted in the last step.) 

In Sec. V, we evaluated the conductance from the P-H 
symmetry violating tunneling term using the Landauer 
formalism. It is instructive to evaluate the conductance 
instead using the Kubo formula, and the commutator 
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approach, as done earlier in this Appendix. In this case, 
the new term in the (a = 0) current operator is: 



6I = t\e l7T{N - 1)/2 ^ + ifj- + h.c 



(A36) 



The Green's function for I , for large imaginary times, 
now picks up a correction: 



5 < Io(r)J (0) >=4A 2 



1 



(2/3) 2 sin 2 (ttt//3) 



(A37) 



We may evaluate the Fourier transform using the same 
methods as above. However, we focus immediately on 
the T = limit, for simplicity: 



SG(u n ) = -4A 2 



00 dre iu} - T 



(2ttt) 2 



(A38) 



This integral requires an ultraviolet cut off, as before. 
Cutting off the integral at |r| > To, gives: 

0\2 poo j 

SG = ^~ ^cosKr). (A39) 

7T Jr„ T 



It is convenient to integrate by parts, giving: 
2A 2 , 2\ 2 uj n f°° dr 



sg 



7r 2 r 



X 2 u n f°° dr . . s fkln . 
» / _ s in( W „r). (A40) 

71 Jr n r 



[We used cos(w„To) ~ 1 in the first term.] Since this 
integral now converges, we may take tq — > 0, assuming 
that |w„|to <C 1. We may then rescale the integral giving: 



sg 



-2A 2A 2 |w. r , 



ds 



7T 2 T0 7T 2 J Q S 



Evaluating the integral gives: 

-2A 2 2A 2 U 

oy 



tt 2 tq 2tt 



sins + O(r ). (A41) 



O(r ). (A42) 



Note that we again obtain a term which is even in u> n 
but singular, involving an absolute value. Analytically 
continuing to real frequency, using Eq. (|A19|) we obtain 
the zero temperature DC conductance from Eq. (|A10|) : 



A 2 2e 2 

TT h 



(A43) 



the same result than is obtained from the Landauer for- 
mula. (Again a factor of e 2 /h was inserted at the last 
step.) 

Alternatively, we could cut off the integral in Eq. 
(IA38|) . by a finite bandwidth, 2D. Then the fermion 
propogator is modified to: 



<Vf(r)Vto(o) >-a^i 



l-e- g H 

2-KT 



giving: 



SG(u n ) = 4 A 2 



(2ttt) 2 



(A44) 



(A45) 



where the integration region is now the entire real line. 
This can be expressed in terms of the exponential inte- 
gral function, Ej(x), giving a term of 0{D), the same 
universal term, — A 2 |w n |/7r, plus terms that vanish when 
l/D — > 0. This confirms that the term cx \u) n \, which 
determines the conductance, is indeed universal. 



APPENDIX B: DEFINITIONS OF T K 



The definition of Tk used in this paper, and defined in 
Eq. (|2.16j) . is the one first introduced by Nozieres^ in 
the original paper on local Fermi Liquid Theory for the 
Kondo model. It corresponds to an arbitrary choice but 
does have the advantage that the coefficient of T 2 /T|- 
in the conductance has the relatively simple value of 7r 2 . 
[See Eq. (|4.20|) .] Another popular definition is the one 
adopted by Wilson, which we refer to as Tjf . This is 
fixed by the requirement that the impurity susceptibility, 
at T > Tk have a particular form: 



Ximp * 

4T 



1 - 



(1/2) In [HT/TD] 



In 2 (T/Tg) 



-O [In 2 [ln(T/T^)] / In 3 (T/T^ 



At low temperatures, 



(gfi B ) 2 w 2 
^ K 



Ximp 



(Bl) 



(B2) 



where the Wilson number, w has the value: 

w = e c+1 / 4 A 3/2 « .4128 (B3) 

and C is Euler's constant (rj .577216). These definitions 
of Tk are related by: 



Tk = (nw/4)Ti 



K ■ 



(B4) 



Rather accurate results have been determined for the 
conductance at all temperatures using numerical renor- 
malization group^l and integrability^, These authors 
generally use a Kondo temperature, T^, related to the 
other ones by: 



Tk = T^/w = (n/A)T K , 



(B5) 



motivated, perhaps, by the fact that the zero tempera- 
ture impurity susceptibility now takes the simple form 
Ximp — » (5A*b) 2 /(4T^). Some experimental papers^ 
on the Kondo effect in embedded quantum dots define 
a Kondo temperature by the condition that the conduc- 
tance take half its zero temperature value at T = T K : 



C(T = T K )/C{T = Q) = 1/2. 



(B6) 
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The theoretical results in [38| indicate that Tjj- w T^l, In 
[32| we wrote the Fermi liquid interaction: 

Hi„ t = -A (V^e) . (B7) 

(where A is not to be confused with the coefficient of the 
particle-hole breaking term in the current paper) and we 
normalized our fermions in an unconventional way so that 
if>AL = v / 27T'0- Thus: 

A =a!v « B8 » 



Numerous other definitions of Tk are also in use in- 
cluding, for example, the frequency scale at which 
lmT(oj,T = 0) is reduced by 1/2 from its maximum at 
zero frequency. 
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